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Can stellar wobble in triple systems mimic a planet? 
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Abstract. The first extrasolar planets have been detected by the measurement of the wobble of the parent star. 
This wobble leads to the periodic modulation of three observables: the radial velocity, the position on the sky and 
the time of arrival of periodic signals. We show that the same wobble, and therefore the same modulation of the 
three observables, can be due to the presence of a more distant binary stellar companion. Thus, the observation of 
the wobble does not, by itself, constitute a proof of a planet detection. In particular, astrometric confirmation of 
a wobble does not necessarily provide a sufficient proof of the existence of a planet candidate detected by radial 
velocity. Additional conditions, which we discuss here, must be fulfilled. We investigate the observed wobble for 
the planet candidates already detected and we find that, for each case, a wobble due to a binary stellar companion 
can be excluded. 

But for apparent Saturn-like planets in wide orbits, there may be an ambiguity in future detections, especially in 
spaceborne astrometric missions. We conclude that, in some cases, a definitive proof for the presence of a planet 
requires further observations such as direct imaging. 
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1. Introduction 

The detection of the first extrasolar planets rests on an 
indirect method, namely the measurement of the reflex 
motion of the parent star. In cases where only the wobble 
is detected, one can ask whether the detection of radial 
velocity (RV) variations 1/ are indeed due to a stellar 
wobble and not to other effects (such as stellar rotation 
or variable stellar activity) 2/ are indeed due to a planet 
and not to other dynamical effects 3/ that the companion 
is indeed a planet (and not for instance a planetary mass 
black hole or strange matter object). In a few cases, the 
planet detection is confirmed by (or was preceded by) the 
detection of a transit of a planet, but the question of the 
planetary explanation of the wobble remains a priori open 
for the other candidates. 

Here we consider the case where the wobble is real 
but due to the perturbation by a distant binary star first 
suggested by Schneider (| 1999(1 . 



2. An approximation: restricted 3 body problem 
for hierachical systems 

Consider a triple hierachical system consisting of a binary 
system (with masses Mi = M2 = M; separation between 
components 2 as) and a third companion orbiting the cen- 



ter of mass of the system in a bigger circular orbit (mass 
M3; radius of the orbit a a)- The perturbation caused by 
the binary system on the orbit of the third companion 
can imitate the perturbation caused by a planet around 
the latter star. In the appendix we have derived the equa- 
tions of motion of such an approximation. The result of 
the perturbation of the binary system is an elliptical pe- 
riodic motion superimposed to the bigger orbit of the third 
star. Studying this perturbation in the plane of the orbit 
of the third star, let the X axis lay from the center of 
mass of the binary system to the target star, and the Y 
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Fig. 1. Orbital elements of the system (not to scale) 
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axis perpendicular to it (see fig ^ , the magnitude of this 
perturbation (see appendix) in both axis is: 



4.5- 



(1) 



Our approximation takes as starting points the sizes 
of the orbits of the binary system {as) and the third star 
{a a) together with the mass of the stars in the binary 
system (M) and gives the amphtude of the perturbation 
{5x\ Sy) which the binary system causes in the movement 
of the third star. 

Stellar wobble and planetary companion produce, in 
principle, the same perturbation: a periodical elliptical 
motion which can be measured either by radial velocity 
or astrometry. 

We have made simulations of such triple systems with 
the code KAPPA. Taking as initial parameters as, o-a 
and M (see figure we can calculate the initial position 
and velocities for the three bodies which are required by 
the code (using lA.ll IA.5I and IA.7|I . Finally, we compare 
the results obtained in the simulation with those expected 
according to lA.6l 

As an example, we are going to choose ub = 
1 AU; aA = 90 AU; M = IMq. In figElwe see the radius 
of the orbit of the third star around the center of mass of 
the system. In a period of three years, it oscillates three 
times (half the period of the binary system, as expected) . 

The radius fits very well toroi 



radius of the orbit ot the third star 
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Fig. 2. Variation of the distance from the third star to the 
center of mass of the binary system. 



3. Imitating a planet 

The perturbation of the binary system can induce a mo- 
tion in the target star imitating a planet. In other words, 
we see a star with a periodic elliptical wobble. We suppose 
that this wobble is caused by the orbit of the star around 
the center of mass of the system star-planet. However, 
this wobble can be in fact caused by a far binary system 
to which our target star is gravitationally linked. In our 



model, this motion will have the amplitude given byQand 
the following period: 



Pl 



1 _ 2^ 

2 ^'""''" VGM 



(2) 



where M is the mass of each star of the binary system, 
as half the semi-major axis of the orbit of the binary sys- 
tem, as described previously. The subscripts * and pl refer 
to the target star and the supposed planet respectively. 
Other parameters are: 



'GM. 



* r>2 
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\ 47r2 * 
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(4) 
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From these expressions, one can derive the mass M 
required for each star of the hidden binary system to lead 
to the observed values for a pseudo-planet: 



M ^ {'i.br''' Mli" Ml" 

\api 



12/5 



(6) 



4. Application to exoplanets detected by radial 
velocity 

One may wonder if the low amplitude wobble detected 
more than 150 stars (for a permanent update, see 
|http://www. obspm.fr/planets) is due to a planet or to a 
more distant binary system. From the point of view of ra- 
dial velocity measurements, a star is considered as single 
if there is no long term drift in its velocity curve. 

In our model, we take a star suppossed to be single 
with a planet companion but in fact it is a star in a triple 
system. The absence of velocity drift imposes a minimum 
value for the distance of an hypothetical companion. 7 = 
GM/a\ being the acceleration of the target star due to 
a companion at a distance a^i, the velocity drift acquired 
over a time span AT is Ay = 7AT = GM/a\^.T. The 
star is single if AF is smaller than the the observational 
limit. Taking from the last years of radial velocity surveys 
AV^ < 10 m/s and AT = 5 yrs, on gets, for M = IMq, 
aA > 300 pc. 

From equations I A. 71 [21 El and O we obtain a relation 
between the mass of the supposed planet and the period 
of the wobble (with M = = 1 solar mass): 



35 AU 



P. 



8/3 



(7) 



a A J V 1 year / 

In figjSlwe represent this relation for a a — 300 AU and 
a A = 50 AU together with the data for most of the known 
extrasolar planets known nowadays. Clearly, none of the 
planets found risk of being a triple system, because the 
effect that a binary system causes is small. 
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MASS vs PERIOD AT A GIVEN a, 



orbital period (years) 



Fig. 3. The perturbation caused by the binary system in 
the target star will imitate a planet whose mass (for a 
given orbital period) lays in the straight lines, which rep- 
resent distances from the target star to the center of mass 
of the binary system of 50 and 300 AU. 

4.1. Planets in binary star systems 

There are presently around 15 planets detected in binary 
star systems (Eggenberger et al. iSOMJ. One may wonder 
to what extent the companion to the target star hosting a 
planet is in fact a binary system inducing a stellar wobble 
imitating the effect of a planet. 

Since the separation a a, the orbital parameters of the 
planet and the mass of the star (given its spectral type) 
are then known; the mass of the hypothetical binary is 
given by equation El In other words, let's suppose that 
there is no planet and that in fact there is a binary system 
perturbing the motion of the target star. We can calculate 
the mass of this system to see if this hypothesis is reliable. 
From the data of Eggenberger et al. 2004, tabled shows 
in each case the value for M derived from equation El 

The values found for M are aberrant (except in the 
case of 7 Ceph); thus the corresponding stellar wobble 
cannot be explained as due to the binary nature of the 
companion star. In 7 Ceph the value of 2 solar masses is 
not aberrant, however, the mass found for the companion 
is 0.4 solar masses (Dvorak et al. I2()(),S|I . Our model is 
consistent with a binary system of two solar masses each 
orbiting at 2.4 AU. That system has not been discovered, 
observational results do not point in that direction; we 
discourage the existence of the triple system and accept 
the presence of a planet. 

5. Application to future astrometric searches 

One may wonder if in the future this problem can become 
serious with future more accurate measurements of stel- 
lar wobbles. Clearly, from fig. O the wobble indiced by 
a distant binary star can merely mimic low mass planets 
on wide orbits. This configuration escapes planet detec- 
tion by radial velocity but is well adapted to astrometric 
detection. 

In fig^we represent the expected discovery space for 
two spatial missions: GAIA (Sozzetti et al. I2003|) and SIM 



Table 1. Experimental data for planets in binary sys- 
tems. Last column is the mass (calculated with 
for each star of an hypothetical binary system which 
would cause in the target star same wobble as the 
planet (see text). Data taken from Eggenberger et al. 
120041 and from the Extrasolar Planets Encyclopaedia 
(http:/ /www.obspm .fr/planets). 
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3 000 
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3.9 
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7 Ceph 
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2.03 
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20 
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10^ 



semimajor axis of the orbit of the binary stellar system in 



AU 

h 



semimajor axis of the orbit of the planet in AU. 
mass of the planet in Jupiter masses. 

mass calculated with El for each star of the hypothetical 
binary system which would imitate the wobble of a planet. 



(Space Interferometry Mission, JPL). The constraints are 
the limiting resolution of 10 and 2 micro arcseconds re- 
spectively and the lifetime of the missions: three years. 
Both will detect planets in regions still unexplored by ra- 
dial velocity methods. 



DISCOVERY SPACE OF DYNAMICAL METHODS 
orbital period (years) 
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Fig. 4. Discovery space for GAIA, SIM and radial velocity 
missions, together with the expected perturbation caused 
by binary systems at 50 and 300 AU. 

In figElwe compare the expected results for SIM and 
PRIMA (Phase-Referenced Imaging and Microarcsecond 
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Astrometry at ESO VLTI). PRIMA has the same resolu- 
tion as GAIA (10 micro arcscconds) but is not constrained 
by the three years lifetime, having access to a low-mass 
long-period region where this stellar wobble effect will be 
more important. 



DISCOVERY SPACE OF DYNAMICAL METHODS 
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Fig. 5. Discovery space for SIM, PRIMA and radial ve- 
locity missions, together with the expected perturbation 
caused by binary systems at 50 and 300 AU. 

The possibility of stellar wobble simulating a planet is 
small, however, for long periods and if the distance from 
the target star to the binary system (a^) is not that im- 
portant (e.g. 50 AU), those effects will have to be taken 
into account. 



6. Conclusion 

By itself a stellar wobble is not a proof that a planet is 
detected. It is necessary to verify that no far binary star 
generates the wobble or to confirm that planet by transit 
or direct imaging observations. For the presently know 
planets, the explanation by a perturbing binary star can 
nevertheless be ruled out. But the sensitivity of GAIA, 
PRIMA and SIM is such that for some regions of the {Mpi, 
Ppi) plane there can be an ambiguity between a true planet 
detected by astrometry and a wobble induced by a binary 
star. 
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Appendix A: 

The perturbation of the trajectory of a body in a triple system by a distant binary system is a classical issue in the 
3-body problem fRov ll975|l . Here we remind the analysis only for self-consistency of our paper. 

Consider a simplified triple hierachical system consisting of a binary system of two equal mass bodies (Mi — M2 = 
M) in a circular orbit with radius as plus a third companion orbiting the center of mass of the binary system in a 
bigger circular orbit with radius a a 3> as- 

Suppose that the motion of the binary system is not perturbed by the third body (let M3 — 0) and let them move 
in the plane z = 0. As we want a circular orbit with a given angular velocity lu, we find that the equation of motion 
of the bodies Mi and M2 is: 



r Ml 



a B {con ujti + sin ujt j 
—as [cosuti + sinuotj 



(A.l) 



2GM 



I GM 
4«| 



The gravitational force per unit of mass in any point "f* = (x, y, z) caused by these two bodies of mass M is: 



~f = -GM 



r2 

\n\-- 



where ~r*i = — Mi i = 1,2 



(A.2) 



To make the problem even simpler, we are going to suppose that the motion of the third star is in the plane defined 
by the orbit of the binary system. In the plane we keep only the cartesian coordinates (x, y) or their polar equivalent 
(r, 9) with the identities for the unitary vectors: 



r = cos t/ 1 + sm j 
9 — — sin 9i + cos 6 j 



then we can separate lA. 21 into the central force (in the direction of the center of mass of the binary system, parallel 
to f) and the angular force (perpendicular to the former, parallel to 6): 



r 1/2 = 

l^l/2|- 
1 



r Ml, 



r cos 9 i + r sin 9 j 



±aB [cos Lot i + siiiLotj 



[r =F OB cos {ujt — 9)]f^ aB sin {cot — 9)9 



' a% =F 2raB cos {ujt — 9) 

15 3 

1 ± 3 e cos (wt - 6*) + y cos^ {tot - 9) - -e^ 



where e i 



as 



< 1 



= -GM 
2GM 



ri 



_r2_ 

k2p 



9 



l + -e^ + -e^ cos2{Lot-9) 



-e^ sm2{Lot-9)9 



(A.3) 



correct to order 2 in e and which is the same expression we can find in equation 14.25 of (Roy I1979|l (under 
our assumption that the masses of the binary system are equal and the third mass is zero). As T*' = {f ^ r9^^ f + 

(^2r9 + r(fj 9; we finally arrive to: 



2r9 + r9 



2GM 



3GM 



1 + -e^ + -e2 008 2(1^^-6*) 
4 4 



sin 2 {loI - 9) 



(A.4) 



As a solution to the system of differential equations (jA.4|l we propose that the third body is making a circular orbit 
of radius ua around the center of mass of the binary system with angular velocity Q. perturbed by another elliptical 
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motion of semiaxes 6^ and 6y both much smaller than a a- So the equation of motion of this body is, in rectangular 
coordinates: 

x{t) — aACOsftt + Sx cos2u!t 

y{t) — a A sin fit ~ dySm2ujt (A. 5) 

Notice: 

1. this is not the most general solution. In fact, it is closer to a perturbation solution: these equations are solution of 
the differential equation system (|A.4|I . 

2. the angular velocity of the perturbation is —2a;. As the two bodies of the binary system have the same mass, every 
half a revolution of the system the third body sees the same configuration of the binary system, in other words, 
the same configuration of the perturbation. 

3. the Keplerian angular velocity goes as a~^-^. As a^i ^ as ^ a; ^ il. 

Taking this into account we arrive to the following equations to order one in the perturbation (that is, taking into 
account that Sx/aA ^ 1; 5y/aA ^ 1; f^/i^ ^ 1;): 



r = \/ + y'^ aA ^ 5x cos 2Lot 9 = arctan ^ k, Vlt 

f = ^ ^x"^ +y^!v -2ujdx sin 2Lut 6 = ^ ft - 2uj ^ cos 2ujt (A.6) 

Introducing this solution ljA.6|l in the system of equations (|A.4|1 we obtain the following identities: 



. 2r ^ 2GAf / 3 , 9 . 

QA^ + 4:UJ Sx cos 2LLit = — 1 + -e + -e cos 2tjt 

V 4 4 

. 2. . „ 3GAf 2 . „ 

iu) Oy sm 2ijjt — — 2 — £ sm 2a;i 

from wich we arrive to: 

^.=4.5 J ^, = 3^ 1^^ = ^ (1 + 0.75 J) (A.7) 

"-A "'A "a \ "a/ 

The mass of the binary system does not affect the amplitude of the perturbation but does affect the period 
(^wobble ~ (2^)^^ ^ M~°-^). The angular velocity fl is bigger than the Keplerian one, that means that for the same 
mass, the system is less bounded (the total energy E = T — U is less negative). 
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